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^ ■ Abstract 
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These notes address the planar gravitational wave solutions of general relativity in empty 

space-time, and analyze the motion of test particles in the gravitational wave field. Next we 

consider related solutions of the Einstein equations for the gravitational field accompanied by 

long-range wave fields of scalar, spinor and vector type, corresponding to massless particles of 

qv | spin s = (0, i, 1). The motion of test masses in the combined gravitational and scalar, spinor 

** — , i or vector wave fields is discussed. 
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1 Planar gravitational waves 

a. Planar wave solutions of the Einstein equations 

Planar gravitational wave solutions of the Einstein equations have been known since a long 
time @-0- In the following I discuss unidirectional solutions of this type, propagating along 
a fixed light-cone direction; thus the fields depend only on one of the light-cone co-ordinates 
(u, v), here taken transverse to the x-y-pl&nc: 

u = ct — z, v = ct + z. (1) 

Such gravitational waves can be described by space-time metrics 

g^dx^dx 1 ' — — dudv — K(u, x, y)du 2 + dx 2 + dy 2 = — c 2 dr 2 , (2) 

or similar solutions with the roles of v and u interchanged. If the space-time is asymptotically 
minkowskian. With the metric (g), the connection co-efficients become 

r v — k y x — -V v — - k r y — -V v — - k c\\ 



All other components vanish. The corresponding Ricmann tensor has non-zero components 

(4) 



' &,xx, tiuyuy — „ 



Ruxuy — tiuyux ~ i^-.xy 

The only non- vanishing component of the Ricci tensor then is 

t*>uu — 7) \-*^ ,xx ~v~ -"- .yy ) = 77 ^trans**- • \& ) 

Here the label trans refers to the transverse (x, y)-plane, with the z-axis representing the 
longitudinal direction. In complex notation 

C = x + iy, (,=x- iy, (6) 

the Einstein equations in vacuo become 

iV = & Ktf = 0. (7) 

The general solution of this equation reads 

K(u,(,{) = f(u;0+f(u;{)=J2 / V {tn{k)e^ ku C + Uk)e lku C). (8) 

Note that the terms with n = 0, 1 correspond to vanishing Riemann tensor: R^, VK \ = 0; 
therefore they represent flat Minkowski space-time in a non-standard choice of co-ordinates. 
For this reason we adopt the convention that e = e\ = 0, which is just a choice of gauge. 

b. Geodesies of planar-wave space-times 

We proceed to solve the geodesic equation in the gravity- wave space-time (0) along the lines of 
ref.g: 

x" + T uX "x v x x = (9) 



Here the overdot denotes a proper-time derivative. The proper-time Hamiltonian satisfies a 
constraint imposed by eq.(|^): 

H = g^x"x v 

(10) 
= — uv — K(u, x,y) u + x + y = —c . 

Because the metric is covariantly constant, the hamiltonian is a constant of motion: 

H = 0. (11) 

This can be checked directly from the geodesic equation (H). Also, as v is a cyclic co-ordinate, 
its conjugate momentum is conserved: 

u = 0, (12) 

with the simple solution u = 7 = constant. Again, this agrees with the geodesic equation, as 

there is no non-vanishing connection component in the u-direction: T^" = 0. 

Only the equations of motion in the x-y-plane depend on the specific wave potential K(u,x,y): 



1 -y 2 

= — K^u 2 = ~—K 



,X1 



2 ' 2 

(13) 

1 9 7 2 

y = ~2 K ,v u = —2 K w 

Eqs. (|10|)-([l3|) specify completely the motion of a test particle, with the conservation of H 
taking the place of the equation for the acceleration in the w-direction: 

7V + -y 2 K(u,x,y) — x 2 + y 2 + c 2 . (14) 

If we now add i 2 to the left- and right-hand side, and remember that 

jv — uv = c t — z , (15) 

we can rewrite the hamiltonian conservation law as 

c 2 i 2 +7 2 i<: = c 2 +r 2 . (16) 

Finally, with v = dr/dt = r/i, the equation can be cast into the form 



dt _ l_j^K/l_ 
l -Tr- V l-v 2 /c 2 (17) 

This equation describes relativistic time-dilation as resulting from two effects: 

(i) the usual special-relativistic time-dilation from the relative motion of observers in the rest- 

and laboratory frame, whose time co-ordinates are r and t, respectively; 

(ii) the gravitational redshift resulting from the non-trivial potential K . 

Now from the conservation of 7 — u = ct — z it follows, that 

7 = d (l - ^) , (18) 

with v z = dz/dt. Eqs. (|l7|), ( p~8J ) can then be solved for 7: 

7 2 1 



K , '- v2 / C 2 



(19) 



{l-v z /cf 



Thus, for a paticle starting at rest at infinity in an asymptotically minkoskian space-time, we 
find 7 = c. At the same time we observe that 

1- v 2 /c 2 

h = K + 7, 7^ (20) 

(1 - v z /c) 

is conserved. Now we recall that in our conventions K is at least quadratic in the transverse 
co-ordinates; hence the components x and y of the transverse acceleration vanish for x = y = 0. 
Furthermore K (u, 0,0) = 0, with the result that the origin of the transverse plane moves at 
constant velocity along the z-axis: 

f_ = l-v z /c vg = 1 - 7 2 /c 2 

c 2 1 + Vz/c C l+7 2 /c 2 ' ^ ' 

In particular, the point at rest in the origin moves along the simple geodesic 

a^(r) = (cr, 0,0,0). (22) 

Taking this geodesic as our reference, the solution for the geodesic motion x^(t) of any other test 
particle at the same time presents a measure for the geodesic deviation between the worldlincs 
of the two particles. 

2 Einstein-scalar waves 

Having discussed the planar gravitational waves (0) in empty space we now turn to discuss 
similar unidirectional wave solutions of the combined system of Einstein gravity and a set 
of massless self-interacting scalar fields. The solutions of the inhomogeneous and non-linear 
Einstein equations, with the energy-momentum tensor that of the right- (or left-) moving scalar 
waves, nevertheless turn out to be a linear superposition of the gravitational field of the scalar 
waves and the free gravitational wave solutions discussed in the first paragraph. 

We introduce a set of massless scalar fields a l (x), i = 1, ..., N, taking values in a manifold 
with the dimensionless metric Gij [a] . In four-dimensional space-time the fields themselves have 
dimension [a] = y/E/l; thus, introducing an appropriate length scale 1//, in the context of 
quantum field theory we could write a % = Vhcfrj 1 , with rf{x) a dimensionless field. 

The starting point of our analysis is given by the gravitational and er-model field equations 

D cov a i + r jk i [o]gi"'dn<T : >d l/ a k = 0, 

8nG - ■ (23) 

R^u = — Gij [a] dpa l d u a j . 

Here the covariant d'Alembertian is defined on scalar fields in the standard fashion 

whilst r^ ■ k [a] denotes the Riemann-Christoffel connection in the target manifold of the scalar 
fields. These equations can be derived straightforwardly from the combined Einstein-cr-model 
action, but we will skip the details of that procedure here. Our aim is to construct simultaneous 
traveling wave solutions of the full set of equations ( |23| ) . Such solutions are actually quite easy 
to find. First, the scalar field equation is solved by taking right-moving fields 

a 1 = a l {u), (24) 



with no dependence on any other co-ordinate. Next we substitute this solution of the scalar 
field into the second equation for the corresponding gravitational field. As before, only the 
Mw-component of this equation survives, reading 



1 , „ 8ttG 
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R uu = --AtransK = — dj [a] d u a % d u o } . (25) 



As this is a linear equation, the general solution consists of a linear superposition of a particular 
solution and the general free gravitational wave of the previous section: 

K(u, C, C) = ^? Gn [a] d u a*d u <jJ(( + /(«, C) + />, ()■ (26) 

Now any specific solution a 1 (u) is a map from the real line into the target manifold of the scalar 
fields. Consider the special case that this curve in the target manifold is a geodesic: 

*£ + Ti *L"L-0. (27) 

du 2 3 " du du ' 

Then the quantity 



' = %M^-> (28) 



da 1 da 3 
du du 

generating translations in u, is constant along this curve: dl/du — 0. Moreover, for Euclidean 
manifolds with non-degenerate metric it is positive definite: I > 0. Observe, that for manifolds 
with compact directions (like spheres) the geodesies may be closed; then the corresponding 
scalar field configurations are periodi. 

The special solution for the accompanying gravitational field now becomes 

K scalar (u,x,y) = — j— (x 2 +y 2 ), (29) 

Cr 

to which an arbitrary free gravitational wave solution can be added. In this special case, upon 
inserting K sca i ar into eqs.fll3|) the transversal equations of motion of a test mass take the 
particularly simple form: 

47TG/7 2 ■• 4^G/ 7 2 
x = — x, y = -. — y. (30) 

Thus the test mass executes a simple harmonic motion in the transverse plane, with frequency 

uo = -1 V4nGI. (31) 

The solutions for the coupled Einstein-scalar field equations discussed here are not the only 
ones of interest. For example, the gravitational waves accompanying expanding domain walls 
in a theory with a spontaneously broken global symmetry can be calculated and have been 
discussed e.g. in |J, |). 

3 Einstein-Dirac waves 

In this section we construct wave-solutions for massless chiral fermions coupled to Einstein grav- 
ity. As before the waves are unidirectional, and both left- and righthanded fermion solutions, 
associated with helicity ±1 quantum states, exist. 

To treat fermions in interaction with gravity, it is necessary to introduce the vierbein and 
spin connection into the formalism. With the local minkowski metric rj — diag(+l, +1, +1, —1), 



the vierbein is a local lorentz vector of 1-forms E a (x) = dx M e°(x) satisfying the symmetric 
product rule 

r, ab E a E b = r, ab e«e b dx>*dx v = g^dx^dx". (32) 

In a convenient local lorentz gauge, the vierbein corresponding to the metric (0) takes the form 



E a = (dx, dy, i {{K - 1) du + dv), ~ ((X + 1) du + d«) 

The inverse vierbein is defined by the differential operator V a = e£d^ such that 

E a V a = dx»8„ 



(33) 



(34) 
(35) 



In components it reads 

V a = (d x , d v , -d u + (K+1) d v ,d u -(K-l) 8 V ) . 
Next we compute the components of the spin connection uj a h = dx^io a b from the identity 

dE a = oj a b A E b . (36) 

With the vierbein (p3j) the spin connection has only one component 



ba 
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-K. x 
V -K, x 







K, 
-K, 



K, x K tX \ 



K, v K y 






o / 



(37) 



In order to construct the dirac operator we introduce a basis for the dirac matrices satisfying 
{j a >J b } — 2ri ab , and define a set spinor generators for the lorentz algebra by a a b = \ [ja,1b\- 
Then the dirac operator is 

1 



7 - J D = 7 a V a --W> 



(38) 



The results we need all depend on the property of the light-cone components of the dirac 
algebra: 

7 « = 7 a e ^ = _ 73 + 7o _ (39) 

This element of the dirac algebra is nilpotent: 



{l u f = o. 



(40) 



The same is true for j v = e®j a = \l u ■ Because of the form of the spin connection (|37j), the 
dirac- algebra valued form u> ah cr ab is itself proportional to 7"; its nilpotency then guarantees 
that the spin-connection term in the covariant derivative (pq) vanishes by itself: 



a be _ 

7 w a a bc 



u be r\ 

7 uj u a bc = 0. 



(41) 



Hence the only vestige of curved space-time left in the dirac operator is the inverse vierbein in 
the contraction of dirac matrices and differential operators: 



7 • D = 7 °V a = Yd» 

( d u - (k - 1) d v 



a\d x +a 2 d v 
V +CT 3 {-d u + (K + 1) d v ) 



-<j x d x - <j 2 d y 
-a 3 (-d u +(K+l)d v ) 

-d u + {K-l) d v 



(42) 
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Here we have introduced the following basis for the dirac algebra: 

*-U IT)-*- 1 -** -°-(o -°i). <«» 

with the CTfe the standard pauli matrices. The zero modes of this operator with the property 
that the energy-momentum tensor only has a non-zero T uu component are flat spinor fields 
ip{u) of the form 

*<•> -*■(*$") = ( -« T) (tO -(-$->)■ (44) 

where xi u ) is a 2-component (pauli) spinor. Indeed, first of all spinors of this type are zero- 
modes of the dirac operator: 

7 • Dip = 0. (45) 



This follows by direct application of the expression ([12|) to the spinor ( [44[ ) , using the nilpotency 
of 7". Moreover, with this property it also follows that the energy-momentum tensor takes the 
form 

V = i ? (7„A, + 1«D») i> = i <5^^7« flu V>- (46) 

To see this, first note that the it-component of the covariant derivative Z? M is the only one that 
does not vanish on ip(u) in general. We then only have to check that in all remaining cases 
with 7^ 7^ 7„ the spinor ip(u) (M4f) gets multiplied by a dirac matrix which can be factorized 
such as to have a right multiplicator of the form 7". Again, as (j u ) 2 = 0, T^ u necessarily is of 
the required form (Hq). 

Finally we remark, that the upper- and lower component of the pauli spinor xi u ) m our 
conventions correspond to a negative and positive helicity state, respectively. Thus we find as 
solutions of the dirac operator in the metric (0) two massless spinor states, corresponding to 
right- moving zero-modes of the dirac operator with helicity ±1, respectively. 

This solution is self-consistent as the only non-zero component of the energy-momentum 
tensor is 

Tuu{u) = -1 [xV] („), (47) 

where the prime denotes a derivative w.r.t. u, and the dagger on \ indicates hermitean conju- 
gation of the 2-component spinor. It is then straightforward to solve the Einstein equation for 
K in the presence of the energy momentum distribution of the spinor field: 

K spmor (u,x,y) = - -^- [xV] («) (x 2 + y 2 ) . (48) 

Again, to this particular solution an abitrary free gravitatonal wave can be added. It should 
be mentioned here, that consistency requires the spinors in the energy momentum tensor (|46|), 
(J47J) to be anti-commuting objects, i.e. if the spinor fields xi u ) are expanded in a fourier series 
of massless matter waves, the co-efficients take values in an infinite-dimensional Grassmann 
algebra. Thus the expression can be given an operational meaning only in the context of 
quantum theory, by performing some averaging procedure. For example, if the spinors form a 
condensate such that the kinetic energy S = — ([x'x']) = constant > 0, then such a condensate 
would generate gravitational waves in which test-masses perform harmonic motion of the type 
(|30|), ( |3f| ) with frequency 

w = -^ V2ttGS. (49) 

c z 



4 Einstein-Maxwell waves 

As the last example we consider coupled Einstein-Maxwell fields. We look for solutions of 
wave-type, using the metric (0). In the absence of masses and charges, the field equations are: 

iV = - ^^ (f^f u x - \g»vF 2 \ , D X F X » = 0. (50) 

With the same metric (0) , we also find the same expressions for the components of the connec- 
tion (|h, and the Riemann and Ricci curvature tensors (0), (pt). Therefore the left-hand side of 
the Einstein eqn. ( p0[ ) is fixed in terms of the potential K(u, y, z). 
As concerns the Maxwell equations, the covariant derivative 

D X F X " = d x F Xfl + V XlJ x F^ + T x /F Xu (51) 

reduces to the first term on the r.h.s., an ordinary four-divergence; this happens because in 
the last term the even connection is contracted with the odd field-strength tensor, whilst the 
middle term contains a trace over an upper and a lower index of the connection, which vanishes 
in our case. 

Thus the Maxwell equation reduces to the same expression as in minkowski space-time, 
and it has the same wave solutions. We consider an elementary wave solution, which in terms 
of the co-ordinate system (||) is described by the vector potential 

A fl = (asinfe(c£-z),0,0), (52) 

with the light-cone components vanishing, and with a a constant transverse vector: a z = 0. 
Of course, arbitrary solution can be constructed from the elementary waves (p2j) by linear 
superposition. With u — ct — z and ui = kc the angular frequency of the wave, the electric and 
magnetic fields are 

Efe(u) = ui&cosku, Bfe(u) = k x a cos/cm. (53) 

As usual for e.m. waves, |E^,(0)| = c|Bj,(0)|, and Ej, • B^, = 0. Indeed, the only non-zero 
components of the full field strength are 

Fm = -Fiu = hat cos ku, i = (x,y), (54) 

all others vanishing. It is now straightforward to compute the stress-energy tensor components 
of the electro-magnetic field, with the result 

T uu = e c 2 k 2 SL 2 cos 2 ku, (55) 

and all other components zero. The Einstein-Maxwell equations then reduce to 

AtransK = ^° k 2 3? COS 2 ku. (56) 

This has the special solution 

K em = ^f- fc 2 a 2 cos 2 ku (x 2 + y 2 ) = ^- E 2 («) CC (57) 

c c 

In view of the linearity of eq. ( |56| ) , the general solution is a superposition of such special solutions 
and arbitrary free gravitational waves of the type (ph: 

K(u,(X) = K em (u,C,C) + f(u,0 + f(u,C)- (58) 



Next we turn to the motion of a test particle with mass m and charge q in the background 
of these gravitational and electro-magnetic fields. These equations are modified to take into 
account the Lorentz force on the test charge: 

5" + r, / /i I/ i A = 2- F%x v . (59) 

With the only non-zero covariant components oiF^ v given by eq.(p4), there are no contravariant 
components in the lightcone direction u. As a result the equation for u is not modified, and we 
again find 

u = 7 = const. (60) 

This also follows, because the electro- magnetic forces do not change the proper-time hamilto- 
nian: 

H = g^yX^x" 

(61) 

• • ft \ ■ 2 i -2 , -2 2 

= = —wo — K(u,y, z)u + y + z = — c , 

except that K(u, x, y) now is given by the modified expression (pq). Therefore v is still a cyclic 
co-ordinate and equation (UJ) for v again follows from the conservation of H: 

jv + j 2 K(u,x,y) = x 2 + y 2 + c 2 . (62) 

As a result we find in this case the same formal expressions for the solution of the equations of 
motion in the time-like and longitudinal directions: 



. dt 1-f^ ,__. 

t= Tr = V 1-vVc 2 ' (63) 

whilst 

1- v 2 /c 2 

h = K +T, 7T2 (64) 

(1 - v z /c) 

is again a constant of motion. In both cases of course K now is the full solution (bq) . 

Manifest changes in the equations of motion are obtained in the transverse directions: 

£ = -— V^K - —kacosku, (65) 

where £ = (x, y) is a transverse vector and V{ is the gradient in the transverse plane. If we 
take for K the special solution ([57]) , we find the conservation law 

We ° k 2 a 2 £ 2 cos 2 ku -\ V '° 2 = h = const. (66) 

c (1 - v z /c) 

Inserting the explicit form of u(t) = 77", eqs.(Bq) then take the form 

I = -^^7 2 /c 2 a 2 cos 2 (7fcr)| - ^fcacos( 7 fcr). (67) 

Equivalently, we can use u instead of r as the independent variable: 



d £ AtteqG , 9 2 o q 

— ^j = 5 — /c a cos (ku) 4 — kacosku. (68) 

ffill^ c 7717 

Clearly, it is useful to decompose £ into components parallel and orthogonal to the electric field 
Efc, which in our choice of electro-magnetic gauge is the same as that of the vector potential a: 

€ = S||+1l> (69) 



with 



It follows that 



£' a o t £ xa (>m\ 



^ £|| 47T£ G ,09 2/, w 9 , 

, ' = = — k a cos (km)4h — feacosfcu, 

au^ <r " 7717 

(71) 

Transforming to the cosine of the double argument, the last equation can be seen to reduce to 
the standard Mathieu equation: 

d^L + 2neoG fe 2&2 + ^ 2fcu) = 

du z c z 

whilst the other equation becomes an inhomogeneous Mathieu equation, with the Lorentz force 
representing the inhomogeneous term: 

d 2 £|| | 2tt£ G ^ 2n ( oj ^ t _ _q_ 

7777 



- 5 — k a (1 + cos 2ku) £n = kacosku. (73) 

du 2 - & " 



Obviously, one may try to find a particular solution to this equation by making an expansion 
in powers of cos ku. The general solution is a superposition of this special one plus the general 
solution of the Mathieu equation (Q) . 

A special case is that of static crossed electric and magnetic fields, obtained in the limit 
k — > 0. Then the eqs.(J72j) and ( |73|) reduce to ordinary homogeneous and inhomogeneous 
harmonic equations: 



^ + ^T- E o^ = 0. 



d% , AirsoG a q 



du 2 X c 4 °^ 
The angular frequency of this harmonic motion is 



(74) 
E . 

777C7 



-^^ Eo = 0.29 x l(T 18 £o (V/m). (75) 

Clearly, the Lorentz force due to the constant electric field produces a constant proper-time 
acceleration of the test charge, but the harmonic gravitational component of the motion is very 
slow for practically realistic electric fields: periods of a year or less require a field strength of 
the order of 10 10 V/m or more. 
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